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Abstract 

Meson corrections on chiral condensate up to next-to-leading order 
in 1/N C expansion at finite density are investigated in NJL model 
with chiral symmetry explicitly breaking. Compared with the mean- 
field results, the chiral phase transition is still of first order, while the 
properties at the critical density of chiral phase transition are found 
to change significantly. 

PACS: 11.30.Rd, 11.10.Wx,12.38.Lg 



1 Introduction 



It is generally believed that QCD undergoes chiral restoration at high tem- 
perature and density. The chiral condensate < qq >, which is regarded as 
the chiral order parameter, behaves quite differently as a function of baryon 
density p as compared to its variation with temperature T . At baryon density 
p = 0, quark condensate < qq > does not change explicitly till T ~ 0.9T C , 
where T c is the critical temperature. However, a model independent study 
[[!]] shows that the quark condensate decreases linearly with p at low density. 
At normal nuclear matter density, p = Po = 0.17 fm~ 3 , < qq > has already 
decreased to about 2/3 of its vacuum value. 

During the last decade, NJL model has been widely used to investigate 
chiral phase transition at finite temperature and density. Most publications 
were based on mean-field approximation, i.e., quark self-energy in the leading 
order of 1/N C expansion and meson in random-phase approximation (RPA). 
At this level, the excited meson modes have no back interaction on quark 
self-energy. A thermal system described by NJL model in the mean-field 
approximation behaves as an effective free fermion gas, i.e., only quarks con- 
tribute to the thermal dynamical potential, while mesons do not play any 
role. This is clearly unphysical, since at least pion should play the dominant 
role at low temperature and density [0. 

Therefore, it is necessary to go beyond the mean-field approximation, 
i.e., considering meson corrections to quark self-energy. In [|[ and |4j], the 
meson effects on chiral condensate at finite temperature have been investi- 
gated. Our work in this paper based on the self-consistent scheme of Q M 
is to investigate the meson corrections on the chiral phase transition at finite 
density with T = 0. 

2 Meson Corrections at T = and /i = 

The two-flavor NJL model is defined through the Lagrangian density, 

£ = Qii^dp - m )^ + G[(^) 2 + (^ 75 tV>) 2 ], (1) 

here G is the effective coupling constant of dimension GeV -2 , and m is the 
current quark mass, assuming isospin degeneracy of the u and d quarks, and 
ip, ip are quark fields with flavor, color and spinor indices suppressed. 
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It is not easy to give the full expressions of quark self-energy and meson 
polarization functions self-consistently. Usually an approximation scheme 
called large N c expansion is adopted. V. Dmitrasinovic et al. proved in |5] 
that quark self-energy E and meson's polarization function IT^ shown in Fig. 
1 are self-consistent to the sub leading order in 1/N C expansion and can keep 
all the chiral relations in the chiral limit. It is clear that the back interaction 
which conserves all the chiral properties is reflected in the contributions from 
the meson propagator to the quark mass, namely in <5E. 

Including current quark mass tuq, the gap equation for quark mass can 
be expressed as 

m = m + Sij + <5£, (2) 

where and are the leading and subleading order of quark self-enery 
in 1/N C expansion, which can be read directly from the Feynman diagrams, 

E^^/A^, (3) 



5S = -8GN c N f m / [—^ — — ( 



(2ir) 8 q 2 — m 2 (p 2 — m 2 ) 2 

3g 2 , , (glT y )\ i 



2\2 



(p 2 — m 2 ) 2 ((p + q) 2 — m 2 ) q 2 — m 2 (p 2 — m 

+ I ^ - 4m2 )], (4) 

(p 2 — m 2 )((p + q) 2 — m 2 ) (p 2 — m 2 ) 2 ((p + q) 2 — m 2 ) 

where gl^g A ^ and ^ are m eson coupling constant in RPA. At this posi- 
tion, it is necessary to point out that to keep the diagrams of meson correc- 
tions to quark self-energy and to meson polarization functions in the next-to- 
leading order in 1/N C expansion properly, we have chosen the internal meson 
propagators in pole approximation. 

Evaluating one quark loop integral, one can get a simple relation between 
the quark condensate and constituent quark mass m 0, 

m-m 

- < « >= -is- < 5 > 
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Meson polarization function II a/ (M = tt, a) can be written as 

U M = u[ PPA) + SU M , (6) 

where Ii^ PA) and 5Hm are pion polarization functions in the leading and 
subleading order of 1/N C expansion, respectively, and 

m M = m < 8 + m ( $ + m$. (7) 

Using the standard way of calculating Feynman diagrams, it is easy to 
write down the RPA polarization functions, 

d A p 1 



u[f A) = 4iN c Nf J 



(2tt) 4 p 2 -m 2 

-2iN c N l( e - f „W, I ^ w — ? ^±- w — ?) (8) 

with t n = and e a = 1. 

After calculating the trace of quark loops, the meson corrections to the 
pion polarization function IT,,-, corresponding to the three Feynman diagrams 
in Fig. |l] can be expressed as: 

ARPA)s 2 



1 



+ 



(p 2 — m 2 )((p + q — k) 2 — m 2 ) ((p + q) 2 — m 2 )((p — k) 2 — m 2 ) 

k 2 (q 2 — eM^m 2 ) 



(p 2 — m 2 )((p + q) 2 — m 2 )((p — k) 2 — m 2 )((p + q — k) 2 — m 2 ) 



. (9) 



1 k 2 (q 2 -e M Am 2 ) 



'((p + ?) 2 — m 2 )((p — A;) 2 — m 2 ) (jo 2 — m 2 ) 2 ((p + g) 2 — m 2 )((p — k) 2 — m 2 ) 
1 2k -q 



(p* — ;//')- t/r — j)i-)((j> -\ </)- — m-)((j> — /,•)' — ;//-) 
fc 2 f/7 2 — f.,4m 2, | 

(10) 



Qo 2 — m 2 ) 2 ((p — k) 2 — m 2 ) (p 2 — m 2 ) 2 ((p + g) 2 — m 2 ) ' 
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with the degeneracy = 3, A CT = 1 and 

A « PA) ) 2 «D 



(27r) 4 q 2 — m 2 (q — k) 2 — m 2 
d 4 p 8mN c Nf(k ■ q — (p 2 — m 2 )) 

(27r) 4 (p 2 — m 2 )((p + q) 2 — m 2 )((p + k) 2 — m 2 



In the similar way, we can get the meson corrections to the sigma po- 
larization functions. However, it should be careful that the upper meson 
propagator for 5U^ in Fig. 1 can only be 7r, but for 5U^ it can be n and a. 

With the above quark self-energy and meson polarization functions, the 
meson mass itim is determined through pole condition 

l-2GU M (k 2 = m 2 M ) = 0, (12) 

and the coupling constant guqq is given by the residue at the pole 

9Mqq = 9k 2 | k 2 =m 2 M . ( 13) 

Another important quantity is the pion decay constant f v which is calcu- 
lated from the vacuum to one-pion axial- vector matrix element. Replacing 
one vertex 275T in n„. by ig 7rqq 'j 5 r y fl T/2, we can get the simple relation 

™lf* = ^0 ^ 

9-nqq 2G 

In the chiral limit, f n satisfies the Goldberger-Treiman relation f % g m q = m 
I- 



3 Meson corrections at finite temperature and 
density 

We now extend the above formulae of SU(2) NJL model beyond mean-field 
approximation to finite temperature and density in the frame of imaginary 
time temperature field theory. 

At finite temperature using the Matsubara formalism and associated finite 
temperature Feyman rules, one arrives at the same integrals shown in the 



5 



last section with the replacement 



/ 



d A p % 




■ a f d 3 p 



(15) 



(2tt) 4 (3 



o (2tt) 3 



for quark integration, with the zero-component of momentum po replaced by 
iu n + fi = (2n + l)inT + /x, and 



for internal meson integration, with the zero-component of momentum go 
replaced by iv n = 2mriT. Here T is temperature, \x the chemical potential, 
and the sums on n run over the Matsubara frequencies u n for quarks and v n 
for mesons. 

While the T = calculation could be performed using covariant, as well 
as non-covariant regularization prescriptions, at T ^ Lorentz invariance 
is always broken and the physical 0(3) regularization restriction \p\ < A 
presents itself most naturally. We take the pragmatic approach of apply- 
ing this at finite temperature, and consider the device of the regularization 
scheme as part of the definition of the model. After the summation over the 
frequencies for quarks and internal mesons, one takes for external mesons the 
extension ia\ — > ko — tum- 

Evaluation of frequency summation in two- and three-loop diagrams is 
not trivial, especially when two frequencies appear in a quark propagator or 
a meson propagator resulted from the energy conservation in Feynman rule. 

4 Numerical results at finite density 

As discussed in the beginning of this paper, the quark condensate behaves 
quite differently at finite density compared with that at finite temperature. 
The meson effects on chiral condensate at finite temperature has been inves- 
tigated in and ||. In this paper, we focus our attention on the meson 
effects at finite density. The second reason to choose numerical study at finite 
density is technical: At T = any frequency sum in imaginary time formu- 
lation is reduced to a step function of chemical potential which simplifies the 
momentum integration greatly. 




(16) 
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As mentioned above, we introduced in our calculations the momentum 
cut-off Ap for quarks, Am for pions and sigmas. We can not determine the 
four parameters for SU(2) NJL model, m Q , G, and the three momentum 
cut-off and Am, since we know only two experimental observables m n = 
139MeV and f w = 93.2MeV and an empirical value of quark condensate 
< q q >= -(250±50MeV) 3 . 

In ||, we have numerically investigated the NJL model parameters by 
fitting pion mass and pion decay constant with an appropriate current quark 
mass rriQ = 5.5MeV. The model parameters in the mean field approximation 
are A F = Q37.7MeV, GA 2 F = 2.16, and - < qq >^ 3 = 2A8.0MeV with quark 
mass m = 330MeV. Going beyond mean-field approximation, we choose two 
series of NJL parameters: A F = 562.3MeV, GA 2 F = 2.31, and - < qq > 1 / 3 = 
223.2Me1/ with A M /m = 1.5 and m = 330MeV, and A F = 51h.7MeV, 
GA 2 F = 2.34, - < qq > 1 / 3 = 218.2MeV with A M /m = 2 and m = 370MeV. 
The degree of meson fluctuations is controled by the ratio A M /m. It was 
found in ||, that for the first group parameters, the meson correction on 
chiral condensate is about 30% , and for the second group parameters, the 
meson correction is about 40%. The big correction in the second case arises 
from the big meson fluctuations characterized by the parameter Am/tu. 

Our results of investigating chiral phase transition by using the above 
parameters will be shown as a function of the scaled density n& = p/po, 
where po — 0.17/m~ 3 is the normal nuclear matter density, and p = | < 
if)ty >= | < 7^7°^ > is the baryon density of the system, which has been 
calculated to the next-to-leading order in 1/N C expansion. The scaled baryon 
density as a function of chemical potential p is shown in Fig. [| the stars 
correspond to the results in mean-field approximation, the pentagrams and 
the solid circles correspond to the results beyond mean-field approximation 
with the first and second group parameters, respectively. 

The jumps in Fig. [| show that chiral phase transition at finite density 
is of first order both in and beyond mean-field approximation. For the first 
group parameter, the lower critical density rib = 0.15 is a little bit smaller 
than that in the mean-field approximation rib = 0.162, and the upper critical 
density n& = 2.1 is a little bit larger than the mean-field value rib = 1-93. 
For the second group parameter, it is found that the lower critical density 
rib — 0.13 is 20% smaller than the mean-field value, and the upper critical 
density rib — 2.79 is 40% larger than that in mean-field approximation. In 
both cases, the meson corrections speed up the chiral phase transition, which 
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is similar to the results at zero density and finite temperature in ||, and 
extend the mixed phase of the first-order chiral transtion, while the degree 
of the change depends on the magnitude of the meson fluctuations. 

We show the quark condensate scaled by its vacuum value < qq > / < 
qq > in Fig. ||. function of the scaled density n& = pj po- It is 

found that the quark condensate decreases to 30% of its vacuum value in 
the mean-field approximation, while it decreases to 10% of the vacuum value 
when meson corrections are included for the first group parameters, and even 
decreases to nearly zero for the second group parameters. This means that 
the meson corrections can help melt the chiral condensate more completely. 

To show the meson effects at low density more clearly, we pick out the 
low density part of Fig. ^. a and enlarge it in Fig. |3|. b. It is found that 
in the mean-field approximation, the chiral condensate and baryon density 
satisfy the linear relation M 

<qq> , SttV 

— = = 1 - 75—5-/0 (17) 

< qq >o f£ m i 

before the phase transition, where the nucleon sigma term = 0.043GeV 
is in good agreement with the experimental value 0.045GeV extrapolated 
from low energy pion-nucleon scatterings. However, when meson corrections 
are considered, this linear relation is more or less broken, depending on the 
degree of meson fluctutions. 



5 Conclusions 

In this paper we extended the self-consistent scheme of SU(2) NJL model be- 
yond mean-field approximation and with chiral symmetry explicitly breaking 
to finite temperature and density, and calculated numericaly the meson cor- 
rections to the quark condensate at finite density. 

The numerical results show that the meson fluctuations can significantly 
change the properties of the system at finite density. Compared with the 
results in mean-field approximation, the meson fluctuations 1) speed up the 
chiral phase transition, 2) extend remarkbaly the mixed phase of the first- 
order chiral transition, 3) help melt the chiral condensate more completely, 
4) break the linear relation between the quark condensate and the baryon 
density. 
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Figure 1: Quark self-energy S and meson polarization function Um in the 
quark and meson propagators. and are the leading and subleading 
contributions to the quark mass. H^ PA ^ and 511^°'^ are the leading and 
subleading order contributions to meson polarization function. The heavy 
solid lines indicate the constituent quark propagator, and the heavy dashed 
lines represent n or a propagator — iD^ PA \q) in RPA approximation. 



10 



a 



3.5 

3 

2.5 

2 - 
1 .5 
1 

0.5 




* * 



0.34 0.36 0.3^ 



0.4 0.42 



Figure 2: The scaled baryon density as a function of chemical potential /j,, 
the stars correspond to the results in mean-field approximation, the pen- 
tagrams and the solid circles correspond to the results beyond mean-field 
approximation with the first and second group parameters, respectively. 
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Figure 3: The scaled quark condensate in the whole density region (a) and 
only before the chiral phase transition (b) as a function of the scaled baryon 
density. The stars correspond to the results in mean-field approximation, 
and the pentagrams and the solid circles correspond to the results beyond 
mean-field approximation with the first and second group of parameters. The 
dashed line is used to represent the linear relation in mean-field case. 
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